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Abstract 

We present an exact calculation for the scattering of light from a single 
sphere made of Faraday-active material, to first order in the external 
magnetic field. We use a recent expression for the T-matrix of a Mie 
scatterer in a magnetic field to compute the Stokes parameters in single 
scattering that describe flux and polarization of the scattered light. 



1 Introduction 

Several reasons exist why one wishes to understand light scattering from a dielec- 
tric sphere made of magneto-active material. Single scattering is the building 
block for multiple scattering. Recently, many experiments such as those re- 
ported by Erbachcr etal. and Rikken etal. ^ have been done with diffuse 
light in a magnetic field. It turns out that the theory using point-like scatterers 
in a magnetic field, as first developed by MacKintosh and John Q, does not 
always enable a quantitative analysis, for the evident reason that experiments 
do not contain "small" scatterers. This paper addresses light scattering from a 
sphere of any size in a homogeneous magnetic field. 

The model of Rayleigh scatterers was used successfully to describe specific 
properties of multiple light scattering in magnetic fields, such as for instance 
Coherent Backscattering, Photonic Hall Effect (PHE) and Photonic Magneto- 
resistance (PMR). The first study of one gyrotropic sphere, due to Ford and 
Werner was applied to the scattering of semiconducting spheres by Dixon 
and Furdyna For the case of magneto-active particles, for which the change 
in the dielectric constant induced by the magnetic field is small, a perturba- 
tional approach is in fact sufficient. Kuzmin |^ showed that the problem of 
scattering by a weakly anisotropic particle of any type of anisotropy can be 
solved to first order in the perturbation. Using a T-matrix formalism, Lacoste 
etal. j^] developed independantly a perturbational approach for the specific 
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case of magneto-optical anisotropy. This was successfully applied to compute 
the diffusion coefficient for magneto-transverse light diffusion . Using the T- 
matrix for a Mie scatterer in a magnetic field we have obtained, we discuss the 
consequences for the Stokes parameters |^ that describe the polarization of the 
scattered hght. 

2 Perturbation theory 

In this paper we set co = 1. In a magnetic field B, the refractive index is 
a tensor of rank two. For the standard Mie problem, its value at position r 
depends on the distance to the center of the sphere |r|, which has a radius a via 
the Heaviside function 8(|r| — a), that equals 1 inside the sphere and outside, 

£(B,r)-I= [(eo-l)I + £F *]e(|r|-a). (1) 

In this expression, I is the identity tensor, Eq = is the value of the nor- 
mal isotropic dielectric constant of the sphere of relative index of refraction m 
(which is allowed to be complex- valued) and ep = 2mVoi?/ci; is a dimension- 
less coupling parameter associated with the amplitude of the Faraday effect (Vb 
being the Verdet constant, B the amplitude of the magnetic field, and uj the 
frequency). We introduced the antisymmetric hermitian tensor — iCijkBk 
(the hat above vectors notes normalized vectors) . The Mie solution depends on 
the dimensionless size parameters x = uja and y = mx. In this paper we re- 
strict ourselves to non-absorbing media so that m and ep are real- valued. Since 
ep ~ 10^* in most experiments, a perturbational approach is valid. 

Upon noting that the Helmholtz equation is formally analogous to a Schrdinger 
equation with potential V(r, a;) = [I — e(B, r)] lj^ and energy w^, the T-operator 
is given by the following Born series, 

T(B, r, u) = V(r, lu) + V(r, ■ Go(c^, p) ■ V(r, + ... (2) 

Here Go(w,p) = l/{iLj'^I — p^Ap) is the free Helmholtz Green's operator in 
gaussian rationalized units for pure dielectric particles, and p = — iV is the 
momentum operator. The tensor of rank two (Ap).y = Sij —piPj/p^ projects 
upon the space transverse to the direction of p. The T-matrix is defined as, 

Tka,kV' =<k,(T|T|k',a' >, (3) 

where | k, a > (respectively | k', cr' >) represents an incident (respectively emer- 
gent) plane wave with direction k and state of hclicity a (respectively k' and 
cr').We will call T" the part of T that is independent of the magnetic field and 
the part of the T-matrix linear in B. We have found the following result 0, 

TL.k'.' =£Fio^< I * I K,.' >- (4) 

where the ^' J ^ (r) are the unperturbed eigenfunctions of the conventional Mie 
problem. This eigenfunction represents the electric field at the point r for an 
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incident plane wave | k, ct > . This eigenfunction is "outgoing" for ^'j^ ^ and 
"ingoing" for ^'j^o-- ^'i- ® resembles the perturbation formula for the Zeeman 
shift in terms of the atomic eigenfunctions, although here it provides a complex- 
valued amplitude in terms of continuum eigenfunctions, rather than a real- 
valued energy shift in terms of bound states. 



3 T matrix for Mie scattering 

In order to separate radial and an angular contribution in Eq. (Q), we used 
a well-known expansion of the Mie eigenfunction ^'j^^. in the basis of vector 
spherical harmonics . We choose the quantification axis z along the magnetic 
field. With this choice, the operator S2, the z-component of a spin one operator, 
can be associated with the tensor — The eigenfunctions of the operator Sz 
form a convenient basis for the problem. The expansion of Eq. in vector 
spherical harmonics leads to a summation over quantum numbers J, J', M and 
M'. The Wigner-Eckart theorem applied to the vector operator S gives the 
selection rules for this case J = J' and M = M' . 

The radial integration can be done using a method developed by Bott etal. 
[ pT| , which gives. 



J,M 



cj Y},M(k) Y}:,,(k') + vj Y}Xk) Y:j:i,(k') 

with the dimensionless parameter: 



(5) 



W = VoBX, 

A is the wavelength in the medium. The meaning of the indices e, m is explained 
in the Appendix. In the limiting case of a perfect dielectric sphere with no 
absorption ( Sm(m) ^ ), the coefficients are given by, 

c^u'^jy ( Aj J{J+1) 



T^^ = ^?f^(----^^ + ^ + ^j]^ (7) 



-d^/ujy f Aj J(J + 1) 
J(J+1) \ y 2/2 

with A J (y) — u'j {y)/uj (y), u,j{y) the Ricatti-Bessel function, and cj and dj 
the Mie amplitude coefficients of the internal field Q . 

Two important symmetry relations must be obeyed by our T-niatrix. The 
first one is parity symmetry and the second one reciprocity. These relations 
can be established generally when the hamiltonian of a given system has the 
required symmetries ( c/. Eq. (15.53) and Eq.(15.59a) P454of Ref. (Q) ). We 
give in the Appendix a less general derivation of these relations for our specific 
problem: 
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T_k(T,-k'(T' (B) — Tk(T,k'cr'(B), 



(8) 



T'-k'-(T',-k-(T( — B) — rkcr.k'cr' (B). (9) 

We emphasize that (t(— k) = — (T(k), i.e. <j indicates in fact heUcity and not 
circular polarization. The helicity is the eigenvalue of the operator S • k. 

3.1 The amplitude matrix 

The amplitude matrix A relates incident and scattered field with respect to an 
arbitrary plane of reference. A common choice is the plane that contains the 
incident and the scattered wave vector, and which is for this reason called the 
scattering plane. We will call the linear base the base made of one vector in this 
plane and one vector perpendicular to it. In this basis, the amplitude matrix 
sufficiently far tor ^ 1^ is simply defined from the T-matrix by. 



Ak.k' = — t;; k' = — a Q • (10) 

is a phase factor that depends on the relative phase of the scattered wave with 
respect to the incident wave and which is defined in Ref. (|^). The complex 
conjugation in Eq. ( [lO| ) is simply due to a different sign convention in Newton 
[|o| . When no magnetic field is applied, the T-matrix of the conventional Mie- 
problem is given by a formula analogous to Eq. (|^) where Cj and Vj have 
been replaced by the Mie coefficients aj and bj, and with M — 1. Because of 
rotational invariance of the scatterer, the final result only depends on cos 0, the 
scalar product of k and k', 6* is the scattering angle. Therefore we get in the 
circular basis (associated with the helicities cr and cr'): 



T^a' = - E 4tTTT + k/4(cos0) + aa'rj^iicose)] . (11) 

j>i •^(•^ + 1) 

Alternatively the T-matrix may be expanded on the basis of the Pauli matrices 

T" = ^ [(^r + ^2*)i + (^r - ^2*)^.] ■ (12) 



In Eq. ([11]), the polynomials ttj^m and tj^m are defined in terms of the Legendre 

J 



polynomials P¥ by 



nj,M{0) = ^/'f (cos 6) , rj^MiO) - ;^^,f (cos 9). (13) 

For M = 1, 7rj_i and tj_i are polynomials of cos 6 of order J— 1 and J respectively 
but not in general for any value of M . When written in the linear basis of 
polarization, Eq. ( pT| ) implies that a Mie scatterer has S3 = Si = as imposed 
by the rotational symmetry. For the backward direction 9 = n, the reciprocity 
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symmetry implies that S3 + S4 = for an arbitrary particle (possibly non- 
spherical) 1^ . We will see that these two properties do not hold anymore when 
a magnetic field is present. 

3.2 General case for when k x k' 7^ 

It remains to express the vector spherical harmonics in Eq. (^, as a function 
of the natural angles of the problem. In Fig. ^ we give a schematic view of 
the geometry. In the presence of a magnetic field, the rotational invariance is 
broken because B is fixed in space. Because our theory treats linear in B, 
can be constructed by considering only three special cases for the direction 
of B. If k and k' are not coUinear, we can decompose the unit vector B in the 
non-orthogonal but complete basis of k, k' and g = k x k'/|k x k'|. This results 
in, 



B=k' 



^1 ^ (B.k)(k.k')-B.k' ^i 

(k ■ k')2 - 1 

^ (B.kO(k-kO-B.k ^i 
(k • k')2 - 1 

+ (B.g)T^^., (14) 
The cases where B is either along k or k' turn out to take the form, 

r^V(B = k) = - (cos 6i)f7 + i?2 (cos 6I)<t'], (15) 

UJ 



T^V(B = k') = -[i?i(cos6')a' + i?2(cos6')a], (16) 



with 



Ri{cos9) = V ^7^— T [C.7^j,i(cos0) -f 2?jT7,i(cos(?)] (17) 

i?2(cos0) - V 4rrru [Djnj^,{cose)+Cjrj^,{cos0)] (18) 

In Ref. we gave an expression for T^^, (B = g) involving a double summation 
over the partial wave number J and the magnetic quantum number M. It 
is actually possible to do the summation over M explicitely, thus simplifying 
considerably the numerical evaluation. Indeed, if one expresses T*^ with respect 
to a z-axis perpendicular to the scattering plane for a given partial wave J, one 
ends up with the following relation between the polynomials 7r,/^j\/ and tj^m, 



7rj,i(cos6l) 2 J2 



.J>M>1 



{J-M)\_ ,^,2 



(J + Af)! 



tj,m(O)'cos(M0) 



+ rj,o(0)2 (19) 
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rj,i(cos0) = 2 J2 



(J + Af)! 



7rj,M(O)'cos(M0) 



+ ^j,o(0)' (20) 



J>M>1 

Upon performing the derivatives of these relations with respect to 9 and com- 
paring to the expression for T^^, (B = g) we find, 



ri,,(B = g) = -{Qm+aa'Q^m (21) 



with 



Qi{9)^-i^Ri{cos9)=isixve-^Ri{cose), 1 = 1,2. (22) 

aO a cos 

We are convinced that a rigorous group symmetry argument exists that relates 

the derivative of T^^,{B = k) with respect to 6 to rJ:^,(B = g). 

3.3 Particular case for when k = k' and k = — k' 

The treatment in section (3.2) becomes degenerate when k and k' are collinear, 
i.e. in forward or backward direction. In these cases, B can still be expressed 
on a basis made of k and of two vectors perpendicular to k. The contribution of 
these last two vectors has the same form as in T^^, (B = g) for 6 = or 6 ^ tt, 
which vanishes. An alternative derivation consists to take the limit — > so 
that i?i = i?2 or 6* ^ TT so that Ri = -R2 in Eqs. (|lHl|). This yields, 

i?i(l) = i?2(l) = V(2J+ 1) {Cj + Vj) (23) 



TT 

J>1 



and 



i?i(-l) - -R2{-1) = -- Y.{-iy+\2J + 1) {Cj Vj) (24) 



TT 

J>1 



This means, 



and 



Ti,k-*-i?i(l), (25) 



LO 



TL-k = *-i?i(-l)- (26) 

UJ 

Both T-matrices contain the tensor $ introduced in Eq. (|^) for the dielectric 
constant of the medium of the sphere. For these two cases, an operator can 
be associated with these T-matrices, which is Sz since we have chosen B along 
the z-axis. For Tj^ j^, the presence of the tensor $ is to be expected since we 
know that the forward scattering amplitude can be interpreted as an effective 
refractive index in a transmission experiment . In the framework of an effec- 
tive medium theory, the real part of Eq. ( p5| ) gives the Faraday effect whereas 
the imaginary part gives the magneto-dichroism (i.e different absorption for 
different circular polarization) of an ensemble of Faraday-active scatterers. 
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4 Magneto-transverse Scattering 



From matrix, we can compute how the magnetic field affects the differen- 
tial scattering cross section (summed over polarization) as a function of the 
scattering angle. Its form can be guessed before doing any calculation at all, 
since it must satisfy mirror-symmetry and the reciprocity relation da/dQ^k — > 
k', B) = d(T / dil.{~'k' — k, — B). A magneto-cross-section proportional to B -k 
or to B • k' is parity forbidden since B is a pseudo- vector. Together with the 
rotational symmetry of the sphere the only possibility is: 

^(k->k',B) =Fo(cos6')-|-det(B,k,k')Fi(cos6i) (27) 

where det(A,B, C) = A ■ (B x C) is the scalar determinant constructed from 
these three vectors. The second term in Eq. (|^ will be called the magneto- 
cross-section. 

The magneto-cross section implies that there may be more photons scattered 
"upwards" than "downwards" , both directions being defined with respect to the 
magneto-transverse vector k x B perpendicular to both incident wave vector and 
magnetic field. An easy calculation yields. 

Act = cr„p - crdow7i = 7^ / dO sin^ 6Fi (cos 9) (28) 

JO 

A non-zero value for Act will be referred to as a Photon Hall Effect (PHE). 

For Rayleigh scatterers, the above theory simplifies dramatically because 
one only needs to consider the first partial wave of J = 1 and the first terms in 
a development in powers of y ( since y ^ I ). From Eqs. (^) and (|^) we find 
that, Ci = — 2j/^/m^(2 -|- m^)^ and Vi = —y^/45m'^. We can keep only Ci and 
drop Vi as a first approximation. Adding all the contributions of Eqs. (^4|) and 
, we find, in the linear base 

^ f tok-k' + itiB ■ (k X k') -itiB ■ k \ 

^^^^'-[ u,B-k' to )■ ^'^^ 

where to = —6inal/uj and ti = —QCiW/lo. This form agrees with the Rayleigh 
point-like scatterer discussed in Ref. p^l . 

A magnetic field breaks the rotational symmetry of the particle. If it is 
contained in the scattering plane, Eq. (^9|) shows that we must have a non-zero 
value for and as opposed to the case when no magnetic field is applied. 
This property still holds for a Mie scatterer, the difference being only present in 
the angular dependence of the elements of the amplitude matrix. A magnetic 
field also violates the standard reciprocity principle as can be seen on Eq. (^. 
This implies that 5*3 -1- 5*4 is non zero for the backward direction 6 = tt. The 
relation 53 -|- ^4 = for the backward direction was derived by Van Hulst, 
but does not apply when a magnetic field is present. In fact the magnetic 
field imposes that 5*3 = ^4 at backscattering. This is readily confirmed by the 
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A A 

k k'= cos 

A 



dk'= d(|) sin a da 



y = magneto-transverse 



Figure 1: Schematic view of the magneto-scattering geometry. GeneraUy, de- 
notes the angle between incident and outgoing wave vectors, (f> is the azimuthal 
angle in the plane of the magnetic field and the y-axis. This latter is by construc- 
tion the magneto-transverse direction defined as the direction perpendicular to 
both the magnetic field and the incident wave vector. Angle a coincides with 
angle 9 in the special but relevant case that the incident vector is normal to the 
magnetic field. 



Rayleigh particle, for which Eq. ( p9|) implies that 5*3 -|- 54 = 2^3 — — 2iiiB ■ k 
for 6 = TT. 

Eq. (H) yields Fi(cos6') ~ VBcos0/k so that Eq. (|2|) gives Act = 0. The 
magneto scattering cross section is shown in Fig. ^ for a Rayleigh scatterer and 
in Fig. 1^ for a Mie scatterer for which a non zero value of Act is seen to survive. 

5 Stokes parameters 

To describe the flux and polarization, a 4 dimensional Stokes vector (/, Q, U, V) 
can be introduced |Q . The general relation between scattered Stokes vector and 
incoming Stokes vector is, 

{I,Q,U,V)=F{I,Q,U,V) (30) 

For a sphere and without a magnetic field, the F-matrix is well known and 
equals 
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4e-06- 



4e-06 




4e-06- 

Figurc 2: Magneto scattering cross section Fi{9) for a Raylcigh scatterer. The 
solid line is a positive correction and the points denote a negative correction. 
The curve has been normalized by the parameter W. No Hall effect is expected 
in this case because the projection onto the y-axis of both corrections cancel. 
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0.02- 




0.02- 

Figurc 3: Magneto scattering cross section Fi{6) for a Mie scatterer of size 
parameter x = 2. The curve has been normahzed by the parameter W. The 
soUd Une is for positive correction and the points are the negative correction. 
Now there is a Hall effect: Aa ^ 0. 
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F° = 





V 



i^ll 









^33 

—F34 






^34 
^33 



(31) 



where 



(32) 



f2 
^12 



Fn = (|5i|2 + |52p)/2 
^^i2-(-|5i|2 + |52|2)/2 

i^33 = (52*^1 + 525r)/2 
^^34 - ti-S*2Si + S2Sl)/2 

Among these four parameters only three are independent since Ff^ 

+ ^"34 • The presence of the many zeros in Eq. (^ij) is a consequence of the 
fact that the ampHtude matrix in Eq. (|o|) is diagonal for one Mie scatterer. It 
is in fact much more general. The form of Eq. (^l|) still holds for an ensemble 
of randomly oriented particles with an internal plane of symmetry (such as 
spheroids for instance) [|l3j. In that case, the averaging is essential to get the 
many zeros in Eq. (^l]). It also holds for a single anisotropic particle in the 
Rayleigh-Gans approximation p4[, Jisf. 

For the Mie case, the anisotropy has two consequences: the F-elements that 
were zero for an isotropic particle may take finite values, and they may depend 
on the azymuthal angle (j). When a magnetic field is applied perpendicular to the 
scattering plane, corrections will appear in the diagonal terms of the amplitude 
matrix. We will use the vector H to denote them. When a magnetic field 
is applied in the scattering plane, the amplitude matrix becomes off diagonal, 
which will fill up the zeros in F'^. We will use the vector G to denote these new 
terms. 

If we call F^ the first-order magnetic correction to the F-matrix one finds. 



F'- = 



H12 
SmG2 



sReGa 
sReG4 
H33 
—H34 



-3mG3 
-5toG4 

H34 

H33 



\ 



(33) 



When B is directed along k, the G terms are given by. 



: {SiRI — 5'2-R2)/2 
: {—S^Rl — >5'2-R2)/2 
: {-SlRl + S2Ri)/2 
= [SlRl + S2Ri)/2 

The general case (forward and backward directions excluded) has 




(34) 



G 



(B-k)(k-k')-B-k' 



(k . k') 



l\2 



G 



B=k 



^ ^ (B.kO(k.kO-B-k ^ 



(k-k 



l\2 



B=k 



(35) 
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and G p.^)^, i s obtained from Gg^j^. by exchanging i?i and i?2 in Eq. (|3^) like in 
Eqs. (na,n6|). Finally, we need, 



H = (B . g;)H3^- , (36) 



with 




-3m(5iQi + 52Q2)/2 
-3m(-5iQi + 52^2)72 
3m(-5iQ2 - 52Qi)/2 
3fJe(5iQ2 - 52Qi)/2 



(37) 



The F-matrix defined in Eq. ( pOD can contain at most 7 independent con- 
stants, resulting from the 8 constants in the amplitude matrix minus an irrele- 
vant phase. Our F-'^-matrix has 12 coefficients (4 for the H vector and 8 for the 
G vector). Therefore 5 relations must exist between these 12 coefficients. These 
relations have not been explicitely derived. 

We can write all the expressions above in a very compact way using the basis 
of the Pauli matrices ||l^ 

4- = irr(^it AO a,) + iTr(AOt a, a,), ^"^^^ 

where Tr is the trace of the matrix, the superscript f denotes Hermite conjuga- 
tion, Oi are Pauli matrices and A^ , A^ are zeroth and first order correction in 
the amplitude matrix defined as the T-matrix from Eq. (p^. 

If the incident light is unpolarizcd, the Stokes vector for the scattered light 
is simply equal to the first column of the F-matrix in Eq. (^3|). For instance 
when B is directed along k, the magnetic field will only affect U — F^i and 
the circular polarization V = F^i which would be zero when no magnetic field 
is applied. We choose to normalize the matrix elements F^j that quantify the 
deviation of the polarization from the isotropic case by the flux Fii without 
magnetic field. In Fig. ^ we plotted these normalized matrix elements for the 
cases where B is directed along k and where B is directed along g. We observe 
that off-diagonal F-elemcnts such as F/2 and Fl^ , are generally more important 
in the angle region of 140 — 170°, and increase with the size parameter. In this 
region, these Stokes parameters seem to be very sensitive to anisotropy as also 
found from studies of Stokes parameters of quartz particles . 

The F-matrix of spherical scatterers in Eq. ( ^ ) contains 8 zeros among its 
16 elements. This property persists for an ensemble of randomly oriented non- 
spherical particles having a plane of symmetry because of the averaging over all 
the orientations. In a magnetic field even spherical scatterers can have a non 
zero value for these 8 elements. Furthermore we have good reasons to believe 
that our theory made for spheres in a magnetic theory should also apply to an 
ensemble of randomly oriented non-spherical particles in a magnetic field, since 
the magnetic field direction is the same for all the particles. 

We have chosen the size distribution [0 and optical parameters of a reported 
experiment ||l8||, when no magnetic field is present , but in which all the matrix 
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F11(B)/F11(B-0) 
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Figure 4: Scattering matrix elements i^g^ (B ^ 0) and Fl^[B)/Fl^{B = 
0) of an ensemble of water droplets as a function of scattering angle for B 
directed along k (top), and Fl^{B)/Fl^{B = 0) and Fl^{B)/Fl^{B = 0) for 
B directed along g (bottom). The refractive index is 1.332, a lognormal size 
distribution has been used with rg// — O.lbfim and cTeff — 0.45 and A = 
632. 8nm. The curve has been displayed for W — 0.1. 



elements of were measured and found to be in good agreement with the 
theoretical evaluation from Eq. (|3l|). For water, the parameter W ~ 2.4 • 10~^ 
for a magnetic field of IT. From Fig. we can therefore expect a modification of 
the order of 2.4- 10~^ in the region near backward scattering for Fg^ (B) /F^^ {B = 
0) when B is directed along k. The magneto-optical effects on polarization are 
very small. Nevertheless they may become significant in multiple scattering, 
which usually tends to depolarize completely the light. 



5.1 Forward and backward directions 



When no magnetic field is present, the situations for 6 — Q ov 9 — it are similar 
because the scattering plane is undefined in both cases. We also have H = 
by Eq. (p6|). The remaining contribution is therefore only determined by the 
G-vector, and the final result reads for 9 = 0, 



I 



Fa^n 



Bk 



y — 3m(z) 
with z = S'i(l)i?i(l). For 9 = n, 



Fo. 



B k 



V 







-3to(z') 






-5Re(z) 







-5Re(z') 






3?e(z) 












3to(z') \ 






(39) 



(40) 



with z' = 5i(-l)i?i(-l). 

The functions -Ri(l) and i?i(— 1) defined in Eqs. ( p7| , p^ are very similar 
to 5'i(l) and Si{—1). Both F-matrices contain only two real- valued indepen- 
dent parameters as do the corresponding T matrices. For unpolarized incident 
light only the Stokes parameter V = Fl^{B) of these matrices is non-zero. 
In Fig. |, all the curves are zero for = and 9 = tt except the one of 
Fl-^{B) / Fl^{B = 0). In other words, unpolarized incident light will produce 
partially circularly polarized light (the degree of circular polarization being pre- 
cisely Fli{B)/Fii{B = 0)) for B directed along k in the forward and backward 
directions. This can be understood from the fact that the effective index that 
one can define from Eq. ( |25| ) suffers from magneto-dichroism [i.e different ab- 
sorption for different circular polarization). 

The modified reciprocity relation in the presence of a magnetic field was 
expressed for the amplitude matrix in Eq. (^. For the F-niatrix it implies 
exactly the different signs in the matrix elements of Eq. ( ^0| ) with respect to 
Eq. m. 
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6 Summary and Outlook 



We have shown that the theory developed for magneto-active Mie scatterers so 
far is consistent with former resuhs concerning predictions of the hght scatter- 
ing by Rayleigh scatterers in a magnetic field. Our perturbative theory provides 
quantitative predictions concerning the Photonic Hall Effect for one single Mie 
sphere, such as the scattering cross section, the dependence on the size param- 
eter or on the index of refraction. 

Using the magneto-correction to the T-matrix we have derived the Stokes 
parameters for the light scattered from a single sphere in a magnetic field. We 
have distinguished two main cases. Either the magnetic field is perpendicular to 
the scattering plane and there will be corrections to the usual non-zero Stokes 
parameters, or when the magnetic field is in the scattering plane, the corrections 
fill up the F-matrix elements which were previously zero. We have discussed 
the particular cases of forward and backward scattering. 

We hope that these results will be useful in comparing them to the situation 
in multiple scattering. Even after many scattering events, we suspect that the 
presence of a magnetic field prevents the Stokes parameters U, V and Q to be 
zero. In single scattering, their order of magnitude is controlled by the parameter 
W — VqBX. In multiple scattering however, this parameter must be replaced 
by fVoBl*, where / is the volume fraction of the scatterers and I* the transport 
mean free path, and fl* ^ 1. We expect to find more significant effects in this 
case. 

We thank Geert Rikken and Joop Hovenier for useful comments. We thank 
the referees for their work, and in particular for mentioning the work of Kuz'min 
etal. 



A Appendix: Derivation of reciprocity and par- 
ity relations 

In the indices for polarization in the T-matrix, the state of helicity a is to 
be referred to the direction of the wave vector immediately close to it. In 
Eq. (^,||), r_ko-,-k'cr' for instance really means T_i(;(T(-k).-k'cr'(-k')- To derive 
these equations, we start from Eq. in which we change both incoming and 
outcoming wave vectors into their opposite: 



TLk.,-k'.'= E (-^) «J.AY},,,,(-k).x.(-k)(-k)Y)^V(-k')-X.'(-k')(-k') 



J,M,X 

(41) 

where aj^\ is a well known coefficient, Xcr(k) is the eigenfunction of the operator 
S • k with eigenvalue (T(k) the helicity. S is a spin one operator acting on three- 
dimensional vectors. The summation is to be performed for X — e,m only, which 
are associated with the two transverse components of the given vector spherical 
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harmonics. Yj^j(k) are weh-defined hnear combinaisons of Yj^j(k),Yj^j_]^(k) 
and Yjlj^j^(k) that obeys Q 

Y)jv,^(k)-k = 0, A = e,m. 

We now use the relations, 

Y},M(-k) = (-l)'^+iY}^^(k) 
Y:^^(-k) = (-l)'^YXM(k) ■ ^ ' 

The eigenfunctions Xo-(k) also change under parity since 

X<T(-k)(-k) -x<T(k)(k). 

Because of this additional minus sign, the parities of the vector spherical har- 
monics are in fact. 



The parity symmetry relation of Eq. (|^) follows from the application of these 
relations into Eq. ( ^ ) . The proof of the reciprocity symmetry relation of Eq. 
is similar, where now the following relations are necessary 

Y):m - (-l)'+''Y^,_^, (44) 

for X — e,o,m and 

X^(k)(k) = X-^(k)(k). 
The change of sign of B is provided by the M factor in Eq. (pd]) as surmised. 
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